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STABILITY AND DRAG , CHARACTERISTICS AT MACH NUMBERS 
FROM 0.8 TO 1.5 OF A FREE-FLIGHT MODEL HAVING 
3-PERCENT-THICK, 60° TRIANGULAR WING AND 
HORIZONTAL TAIL SURFACES 
By Rowe Chapman, Jr., and Harvey A. Wallskog 

SUMMARY 

A flight test of a rocket-propelled model having 6o° triangular 
wing and horizontal tail surfaces 5 percent thick was made and data were 
obtained on the drag and stability characteristics. Analysis of the 
data eind comparison with other data on triangular -wing configurations 
indicated that no severe penalty in zero-lift drag was experienced 
because of the- addition. of _a h.o_r_iz^ntal tail and a single vertical tail 
to a low dreig configuration. 

The trimmed drag for a lift coefficient of 0.125 at a Mach number 
of l.li- is 57 percent greater than the zero-lift drag for that Mach 
n\miber . 

The forcing of longitudinal motions by lateral oscillations can be 
a factor in causing a portion of the low amplitude. oscillations in the 
longitudinal plane. 


INTRODUCTION 


A rocket-propelled model having a 5-percent-thick, 6o° triangular 
wing and a 60 ° triangular horizontal tail was flight tested to obtain the 
stability and drag under lifting conditions. The design of the model 
used in this investigation paralleled that of a general research vehicle, 
used for zero-lift drag investigations, on which low drag coefficients 
have been obtained (ref. l) . The primary external difference between 
the present model and the general research model was the presence of a 
single vertical tail and a horizontal tail; hence, the present model 
more nearly represented a possible low-drag airplane configuration. 

A pxxlsing system for disturbing the model in pitch by movement of 
the horizontal tail at a preset frequency was Incorporated for the 
purpose of obtaining the variation of drag with lift over the Mach 
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nvmiber range. Difficulty was experienced with the unconventional 
pulsing system. The final result during flight test was one long pulse 
after which the horizontal tail returned and remained at zero incidence. 
This equipment failure limited the data obtained from the flight test. 

An analysis of the flight time history was made to obtain some of 
the basic aerodynamic coefficients that determine the performance eind 
stability of the configuration. A portion of the flight time history 
evidenced coupled lateral and longitudinal motions; this motion is 
discussed in detail in regard to the frequencies present in the flight 
record. 

The model was flown at the Langley Pilotless Aircraft Research 
Station at Wallops Island, Va. 


SYMBOLS 


Cn 


cc 


^L 

Cd 


Cy 


an 


w 

g 


normal-force coefficient. 

On 

w 


g 

qS 

chord-force coefficient. 


w 


g 

qS 

lift coefficient, Cjg cos 

a - 

Cq sin a 


drag coefficient, Cq cos a + Cjg sin a 
pitching-moment coefficient 

side-force coefficient, — — 

g qS 

yawing -moment coefficient 

normal acceleration as obtained from accelerometer, ft/ sec^ 

longitudinal acceleration as obtained from accelerometer, 
ft/sec^ 

transverse acceleration as obtained from accelerometer, ft/ sec^ 

model weight, lb 

acceleration of gravity, ft/sec^ 
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p 

Po 

V 

1 

M 

1 

5 
R 
c 

X 

A 

I 

lY 

iz 

Ix 

e 

a 

3 

0 

6 
P 


free-stream static pressure, Ib/sq ft 
standard sea-level static pressure, Ib/sq ft 
velocity, ft/sec 

dynamic pressxire, i Ib/sq ft 

Mach number 

specific heat ratio (l.4o) 

wing area (including area within the fuselage), sq ft 
Reynolds number, based on wing mean aerodynamic chord 
wing mean aerodynamic chord, ft 

longitudinal distance along body axis from nose, in. 
radius of equivalent body of revolution, in. 
cross-sectional area in plane normal to body axis, in.^ 
body length, in. 

p 

moment of inertia about Y-axis, slug-ft 

p 

moment of inertia about Z-axis, slug-ft 
moment of inertia about X-axis, sltig-ft^ 

inclination of principal axis to longitudinal axis of model, 
deg 

angle of attack, deg 
angle of sideslip, deg 

angle of attack of principal longitudinal axis of inertia, 

T) = a - e , deg 

angle of pitch, radians 

incidence of horizontal tail, deg 

period of oscillation, sec 
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time , sec 


a 


rate of change of angle of attack, — i— — , radians/sec 

57.3 dt 

rate of change of angle of pitch, d.0 /dt 


Cniq - 


dC, 


m 


di£ 

2V 




dC 


m 


d 

2V 


Cn * static-directional-stability derivative as obtained from the 

P 


formula = 


qSbp2 


T 1/2 time to damp to one-half an^jlitude, sec 

Subscripts: 

T trimmed, or mean value 

The symbols a and 3 used as subscripts indicate the derivative 
of the quantity with respect to the subscripts. 


MODEL AND INSTRUMEINTATION 
Model Description 


The general arrangement of the model is shown in figure l(a) , a 
three-view drawing of the model. Figvire l(b) is an area-distribution 
plot for the model, showing the geometric relationship between the model 
and its equivalent body of revolution. Photographs of the model are 
presented as figure 2. Pertinent geometric and mass characteristics of 
the model are given in table I. 

Fuselage construction was of laminated wood with metal rings from 
station 25 to station 120.5 and housed an ABL Deacon sustainer rocket 
motor, the principal power plant. A spvin metal nose and tail conqileted 
the fuselage. The space inside the nose cone was utilized to house the 
telemetry. Batteries, hydraulic accimixilator, and accelerometers were 
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located in compartments routed into the wooden fuselage and accessible 
through metal hatches. The hydravilic actuator for the tail svirface was 
located inside the hollow rear portion of the vertical tail. 

The wing was of composite wood and metal construction, aluminum 
overlays, and trailing-edge inserts. The horizontal tail sirrface was 
geometrically similar to the wing and was machined of solid aluminum. 
The vertical tail was a constant thickness of 7/8 inch at the juncture 
of the straight trailing edge and the contoured forward portion. 

The pulsing system for the horizontal tail was designed on the 
principle of an unbalanced force existing between a spring and the 
hydraulic fluid acting on a piston. Since no measirrements pertaining 
to the operation of the pulsing system were telemetered, reasons for 
improper operation during flight could not be ascertained. 


Instrimientation 

The model contained a standard NACA eight-channel telemeter which 
transmitted continuous flight measurements of angle of attack, normal 
acceleration at the center of gravity, normal acceleration at the nose, 
transverse acceleration, longitudinal acceleration, horizontal tail 
position, free-stream total pressure, and a reference static pressure. 

Model position in space was determined from an SCR 584 radar 
tracking \mit and model velocity was obtained by use of a CW Doppler veloc 
imeter unit. Atmospheric data were obtained from a radiosonde released 
just prior to model flight.* 


TEST AND ANALYSIS 
Test 


The model was launched at £in elevation angle of approximately 6o° 
utilizing a mobile launcher. Figure 2 (b) is a photograph of the model on 
the launcher and shows the 5- 0-inch high-performance air-to-ground (HPAG) 
rocket motor used for a booster. The booster propelled the combination 
to a Mach number of approximately 0.3 at which time the sustainer rocket 
motor fired, separating the model and booster. The model attained a 
maximum Mach number of I.5I At sustainer burnout. 

The velocity of the model from radar was corrected for the curved 
flight path and this in conj\mction with the atmospheric data from the 
position radeir plots was used to compute the Mach number. The static 
pressure from radiosonde data and radar position plots was used for al ^ 
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data reduction. A check on the Mach number by use of the total press\xre 
and static pressvire was made and a maximum disagreement in Mach number 
of 0.05 was obtained from peak Mach number to a Mach number of l.k2. 

Data corrections were applied in accordance with the accepted procedxires 
as are given in references 2 and 5* 

Figiire 5 gives the variation of the test conditions with Mach num- 
ber for the power-off portions of the flight. In figure 5 (a), the 
Reynolds number is based on the mean aerodynamic chord of the wing. 


Analysis 

The analyses are based on the assvimptlon of constant coefficients 
in the differential equations of motion which in turn are for small 
disturbances from trim conditions. 

Oscillations of the normal acceleration and angle of attack were 
analyzed by use of the method presented in reference k. Oscillations 
of the lateral accelerometer were analyzed by the use of the formula 
presented in reference 5* 

The periodic motion recorded by the normal accelerometer after the 
horizontal tail retvirned to the zero incidence position was not easily 
identified.- A repetitive pattern on the accelerometer trace indicated 
that a steady-state type of motion was occurring. Since the lateral 
motion was of sinusoidal character, any forcing of the motion in the 
normal plane should also have sinsuoidal characteristics. A harmonic 
analysis of the normal accelerometer trace was made in order to separate 
and identify the frequencies present on the record. The nximerical 
method, Rvinge's schedule, presented in reference 6 was used to accomplish 
the harmonic analysis. 


Accuracy 

The telemetered data are believed to be accurate to within ±1 percent 
of the fvill-scale range of the respective instniments. Converted to 
coefficient form, the probable errors in the basic coefficients and angle 
of attack are as follows: 


M 

^N 

^C 

Ax, deg 

0.8 

0.025 

0.002 

0.15 

1.5 

.008 

.0007 

• 15 
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A portion of the errors presented in the aforementioned table are 
introduced by possible errors in dynamic pressure. The error in dynamic 
pressure is estimated to be less than 0.28 Ib/sq in. at subsonic speeds 
and less than O.5O Ib/sq in. at supersonic speeds. The possible error 
in Mach number is of the order of ± 0 . 02 . This error can in paxt explain 
the discrepancy in Mach number between that computed from total pressure 
and that obtained from radar. 


RESULTS AND DISCUSSION 
Performance 


Figure presents the variation of the horizontal tail incidence 
with Mach number. The solid line is for power-off flight and the broken 
line is for power-on flight; below M = I.06 the tail incidence remained 
at zero for the power-on condition. Analysis of the data was made from 
various portions of the record where the model maneuvered. 

Longitudinal trim .- Figure 5 shows the trim characteristics of the 
model. The solid lines are the Cl and a trim for the power-off 
portion of the flight. The broken line is the trim angle of attack for 
the power-on condition, and, at Mach nijmbers less than 1 . 1 , it represents ' 
the trim attitude necessary to counteract the pitching moment caused by 
the rocket motor. The power-on curve is shown for the purpose of 
relating the attitude of the model to the lateral stability during the 
accelerating portions of the flight. 

The changes in trim angle of attack for the power-off condition are 
a primary result of the changing incidence of the horizontal tail. The 
changes in trim lift coefficient as contrasted to the changes in trim 
angle -of -attack are additionally affected by a changing lift-curve slope. 
The transonic region was traversed with no abrupt changes in trim for 
the low angles of attack of this test. 

Lift . - Figure 6 is a plot showing the variation of Cl and Cjj 

with a during the two maneuvers caused by the moving horizontal tail. 

The Mach ntimber and horizontal -tail position corresponding to the lift 
and angle of attack are shown in the figure. 

Figure 7 is a plot of the slopes of the Cl against a curves as 
a function of Mach number. Additional points of Cj;^ were obtained 

from low-range a variations at the subsonic Mach numbers anti are also 
shown in figure 7- 
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For comparison a ciirve of Cj^ for an airplane configuration 

having a 60° triangvilar wing but no horizontal tail is shown in figure 7. 
This ciurve is from reference 7 and is for a model with a 6.5-percent- 
thick wing, 6=0^ and power-off conditions. An additional comparison 
is provided by the curve labeled reference 8, which is for the wing plus 
wing-body interference. This c\irve is for configuration A from refer- 
ence 8 and these data agree very well with the low-lift data from the 
present test. The supersonic Cj^ seems to agree in level with the 

data from reference 7; however, the transonic level is somewhat higher. 

Drag .- Presented in figure 8 are the power-off drag coefficients 
plotted as a fxmction of lift coefficient for M = I.52. The lift 
coefficient and tail-incidence values correspond to those shown in 
figure 6(b) and additional points for M = 1.49 are included. 

A single value of the induced drag parameter dCp/dCL^ =0.51 
at M = 1.52 was calculated by utilizing the cxirve presented in figure 8. 
This value of induced drag corresponds to the trimmed value since the 
tail deflection 6 was changing slowly. Untrimmed induced drag 
measurements were presented in reference 7 for a delta-wing configuration 
where the eleven controls were changed abruptly from 6=0° to 6 = -9° 
and in reference 9 for a delta-wing model which had an unswept horizontal 
tail that was pulsed abruptly from 6 = -l.l6° to 6 = -5.45°. From 
these data estimates of the trimmed dCp/dCL^ were obtained. 

The \mtrimmed value of dCp/dCL^ for the tailless delta-wing model 
of reference 7 (6.5-percent-thick wing) varied from approximately 
0-35 (b = 0°) to 0.40 (& = -9°) and the estimated trimmed value was 
about 0.6. For the delta-wing model with an unswept tail (ref. 9), the 
untrimmed dCp/dCL^ varied from about 0.2 (d = -l,.l6°) to O.3 (8 = -5.45°) 
and the estimated trimmed value was approximately 0.4. Hence-, it appears 
that, althoiagh the induced drag obtained from the present model was high, 
it con^jares favorably with other tests, since it is lower than the 
corresponding value estimated for a tailless model and higher than that 
obtained from a model which had an unswept tail that was relatively 
larger and placed farther rearward. 

The maximum L/D obtained from the present test was 4.4 at a 
value of Cl of 0.12 compared with 7*0 reported in reference 1 for a 
model with the same wing-fuselage combination but with a different 
empennage. The low value of (L/o)j„a^ from the present test was a 

result of the high trimmed induced drag, whereas the reference data 
utilized untrimmed values of dCp/dCL^. 
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Illustrated in figure 9 are the measured drag coefficients plotted 
against Mach number. Corresponding trim lift coefficients are shown in 
figure 5(a) . The point on figxire 9 labeled zero-lift drag was obtained 
by extrapolating the drag polar of figure 8 to Cl = 0. At Mach numbers 

below about 1.5 the measured drag data presented correspond very nearly 
to zero-lift drag. The dashed line in figure 9 labeled reference 1 is 
the zero-lift drag of the general research vehicle which was the same 
configuration as that of the present test except for the empennage. The 
solid line in figure 9 is the estimated zero-lift drag for the present 
test and was obtained by drawing the curve through the zero-lift drag 
point (from the drag polar) and peirallel to the drag curve of the refer- 
ence model. From figure 9 it is apparent that the zero-lift drag 
penalty caused by the change in the empennage was about 0.0020. Fin 
drag for the reference model was experimentally determined to be about 
0.0012 at supersonic speeds. The increments in drag from the estimated 
zero-lift drag and the measured drag points are indicative of the drag 
due to lift for 6 = -8° of the present model. The trimmed drag and 
lift coefficients which resulted from a tail deflection of 6 = -8° at 
a Mach number of l.ll- were 0.0282 and 0.125, respectively. This drag 
increment corresponds to 57 percent increase in drag to obtain a lift 
coefficient of 0 . 125 . 


STABILITY 

Longitudinal 


Three primary motions occurred during the flight that provided data 
^imenable to the linear analysis technique. When the sustainer motor 
fired, a transient oscillation was induced because of the pitching moment 
incurred from the thrust line not passing through the center of gravity 
of the model. From this oscillation the period of the model was measured 
and adjusted to the power-off Inertia of the model. This period is 
plotted in figure 10 and is labeled power on. When the motion of the 
horizontal tail occurred, a transient oscillation followed and the period 
of this oscillation is plotted in figure 10 and is labeled as the power- 
off period. 

Static . - The values of reduced from the period curve shown 

in figure 10, are plotted as a function of Mach number in figure 11. 

Data from reference 7 (for a tailless configuration having a different 
fuselage) are presented for coir^iarison. These data from reference 7 for 
center-of -gravity locations of 0.20C and 0 . 25 c bracket the cxorve obtained 
from the present test which is for a center-of -gravity location of 0 . 257 c. 
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Dynamic Two data points of 


Cmq + 


were obtained from 


analysis of the dan^jing of the transient oscillation. These data points 
and a comparison with data from reference 7 are shown in figure 12 . The 
con5>arison shown in figirre 12 indicates that the horizontal tail was 
effective for increasing the damping of a triangular -wing configuration 
at M = 1.5. 


Lateral 

For the power-off condition and the trim attitude shown in figvire 5 
the model experienced dynamically unstable, neutrally stable, and stable^ 
conditions in traversing the Mach nvimber range. The model oscillated in 
yaw from M = O.75 to M = I.3 where the trim attitude of the model was 
changed by the horizontal— tail motion and the model bec6une dynamically 
stable. The data obtained from the model oscillations in yaw at the 
power-off conditions are presented in figures I5 to 16. 

The measured period as a function of Mach number, presented in 
figure 15, shows an abrupt shift at M = O.96. This change in period 
is more evident and distinctly present on a plot of period against flight 
time. 


Static .- The static-directional-stability derivative 


obtained 


from the measured period, is shown in figure 14 as a function of Mach 
number. This curve indicates that a rearward shift of the lateral 

aerodynamic center occurs in going from subsonic to supersonic Mach 
numbers. Also indicated by the ^np* curve is a stability increase 

between M = 1.0 and M = 1 . 15 ; this increase is in a region which 
corresponds with that in which an increase in lift-curve slope for sur- 
faces normally occiirs. The Cnp for the wing-fuselage combination, 

obtained from reference 10 for zero a, is shown in figure lit-. The 
abrupt change in period and reflected change in C^p* which occurred 

at M = 0.96 corresponds to a changing lateral trim at the same Mach 
number . 


A computation was made to determine the difference in Cnp* 


obtained from a one-degree -of -freedom analysis and a reduction utilizing 
estimated derivatives by the three-degree-of -freedom assun5)tion. The 
results indicate that Cnp* is approximately 8.1|- percent lower than Cnp. 

References 11 , 12 , and I3 were used for estimating the derivatives and 
for computing the lateral motions. 
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Dynamic For the power-off conditions the lateral dynamic stability 
as a fvmction of Mach number, given in figure I 5 , shows a region of 
neutral dynamic stability from M = 0.75 to M = 0.88, a region of 
dynamic stability from M = 0.88 to M = 1.0, and a region of dynamic 
instability from M = 1.0 to M = 1 . 3 . Between M = I .3 and M = I. 5 , 
the model was trimmed to a higher angle of attack and was apparently 
dynamically stable. 


For the power-on conditions and the model trim as shown in figure 5 , 
no dynamic instability was apparent as the model traversed the Mach 
nvunber range in accelerated flight. 

Apparently two factors were of primary influence in the dynamic 
stability for the power-on conditions. The inclination of the principal 
axis was favorably changed by the additional mass of the rocket grain. 

The power-on trim attitude of the model was more positive below M = O. 90 . 
Other factors believed to be minor are jet effects, small center-of- 
gravlty change, effect of angle of attack on the derivatives, and so 
forth. A plot of T] as a function of Mach niimber for the power-on and 
power-off conditions is shown in figiare I 6 . 


Lateral trim .- The trimmed side-force coefficient of the model is 
shown in figure I 7 . A significant factor of the lateral trim character- 
istics is that a change in lateral trim occvmred at M = 1.33 when the 
inclination of the horizontal tail was changed. This trim change, which 
corresponds to a sideslip angle of approximately 1°, is not large in 


magnitude but introduced a transient lateral disturbance for rapid motions 
of the horizontal tail. Another change in Cy trim occurs between 
M = 0.88 and M = O. 98 . This transonic lateral trim change is estimated 
to be of the order of O.L° sideslip. 


Coupled Motion 

Figure I 8 presents a typical time history of the model flight for 
the low lift attitude, power-off condition. The sustained motion in the 
transverse direction can be explained by a dynamic instability for the 
low lift attitude. The motion in the normal direction as shown by the 
normal accelerometer and the angle -of -attack measurements is of the 
order of sustained angle -of -attack disturbance. 

In order to determine the cause of the disturbance in the normal 
plane, the angle-of -attack motion was broken down into its two harmonic 
conponents. Figure 19 shows the experimental data points taken from the 
flight time history at a Mach number of 0.92. The two sinusoidal compo- 
nents, curves A and B, can be added to give the conpoxinded motion which 
agrees very well with the experimental data. Curve A has the same 
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frequency as the lateral motion but there is a phase shift between this 
motion and. the motion indicated by the transverse accelerometer. Curve B 
has. a frequency twice that of the lateral motion and is also shifted by a 
phase angle. 

The occurrence of motions in the normal plane that have frequencies 
twice that of the lateral motion can be satisfactorily ejfjlained by 
either the derivative Cjjj^ or by terms present in the equations of 

motion before linearization, such as the pi^ term in the equation for 
pitching motion. However, the frequency of the motion labeled curve A 
is different from the nattiral pitch frequency of the model (P = O. 51 ) . 

The time history shows the horizontal tail to be oscillating, but only 
within 2 percent of the full-scale range of 8 °. The correspondance of 
one pitch-plane frequency. with the natural frequency in yaw remains 
unexplained and what, if any, part the elastic (spring and mass) control 
system plays in causing the oscillation is unknown. 


CONCLUSIONS 


A large-scale rocket-propelled model was tested in free flight at 
Mach nmbers up to 1 . 5 . The model which represented a possible low-drag 
airplane configuration had a 5 -percent-thick, 6o° triangular wing and a 
geometrically similar horizontal tail located on top of a sweptback 
vertical stabilizer. The horizontal tail was located approximately at 
50 percent mean aerodynamic chord above the wing chord plane anH hari a 
tail length of 125 percent mean aerodynamic chord. From analysis of 
the flight test data and comparison with data from similar configurations, 
the following conclusions are presented; 

1. No severe penalty in zero-lift drag was experienced because of 
the addition of a single vertical and horizontal tail to a low drag 
configuration. 

2. The trimmed induced drag parameter of O .51 obtained from the 
present test at a Mach number of 1.52 and values of lift coefficient up 
to 0.15 was high in comparison with the untrimmed values obtained from 
tests of similar configurations with and without tails but compared 
favorably with estimated trimmed values. 

5 . The trimmed drag and lift coefficients which resulted from a 
tail deflection of - 8 ° at a Mach nvmiber of 1 .^ 4 - were 0.0282 and 0.125, 
respectively. 
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1|-. The forcing of longitudinal motions by lateral oscillations can 
be a factor in causing low ainplitude oscillations in the longitudinal 
plane . 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va., July 2, 195^. 
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TABLE I 

PHYSICAL CHARACTERISTICS OF THE MODEL 
Geometric 

Wing: 

Area (included), sq ft 15. 9 

Span, ft 5.75 

Aspect ratio 2.5I 

Mean aerodynamic chord, ft 5.5O 

Sweepback of leading edge, deg 60 

Dihedral q 

Taper ratio 0 

Airfoil section NACA 65AOO5 

Vertical tall: 

Area (included), sq ft . I1..65 

Height (from center line), ft I.67 

Aspect ratio O.6O 

Sweepback of leading edge, deg 66.7 

Taper ratio 0.28 

Horizontal tail: 

Area, sq ft 1.59 

Span, ft 1.92 

Aspect ratio 2.5I 

Sweepback of leading edge ^ deg 60 

Taper ratio ^ q 

Airfoil section NACA 65AOO5 

Mass and Inertia 

Fully Grain 

loaded expended 

Moment of inertia: 

About X-axis, slug-ft^ 5.5 5.57 

About Y-axis, slug-ft^ 85.9 68.2 

About Z-axis, slug-ft^ 85.7 68 

Principal axis inclination, deg 1.82 

Weight, lb ^ 451.0 555.6 

Center of gravity, percent C 21.7 25.6 
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Figure 1.- Geometric ' characteristics of model. 
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(b) Area distribution and equivalent body characteristics. 
Figure 1.- Concluded. 
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Figure 2.- Photographs of model. 
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(b) Model on laiincher. 
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Figure 2.- Concluded. 
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Figure 4.- Horizontal-tail position. 
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(a) Trim lift coefficient. 
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(b) Trim angle of attack. 

Figure 5.- Model trim characteristics. 
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Figure 10.- Period of normal oscillation. 



Figure 11.- Static longitudinal stability. 
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Figure 12.- Dynamic longitudinal stability. 
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Figure l4.- Lateral static stability. Power off. 
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Damping factor (-I/T1/2) 











Compounded motion 


Experimental data 
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